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Rings are a model system for studying phase coherence in one dimension. Superconducting rings 
have states with uniform phase windings that are integer multiples of 2n called fluxoid states. 
When the energy difference between these fluxoid states is of order the temperature so that phase 
slips are energetically accessible, several states contribute to the ring's magnetic response to a flux 
threading the ring in thermal equilibrium and cause a suppression or downturn in the ring's magnetic 
susceptibility as a function of temperature. We review the theoretical framework for superconducting 
fluctuations in rings including a model developed by Koshnick^ which includes only fluctuations in 
the ring's phase winding number called fluxoid fluctuations and a complete model by von Oppen 
and Riedel^ that includes all thermal fluctuations in the Ginzburg-Landau framework. We show 
that for sufficiently narrow and dirty rings the two models predict a similar susceptibility response 
with a slightly shifted T c indicating that fluxoid fluctuations are dominant. Finally we present 
magnetic susceptibility data for rings with different physical parameters which demonstrate the 
applicability of our models. The susceptibility data spans a region in temperature where the ring 
transitions from a hysteretic to a non hysteretic response to a periodic applied magnetic field. The 
magnetic susceptibility data, taken where transitions between fluxoid states are slow compared to 
the measurement time scale and the ring response was hysteretic, decreases linearly with increasing 
temperature resembling a mean field response with no fluctuations. At higher temperatures where 
fluctuations begin to play a larger role a crossover occurs and the non-hysteretic data shows a fluxoid 
fluctuation induced suppression of diamagnetism below the mean field response that agrees well with 
the models. 

PACS numbers: 73.23.Ra, 74.40.-n, 73.23.-b, 74.25.Ha, 74.25.Dw, 74.62.-c 



I. INTRODUCTION 

Fluctuations play an important role in the supercon- 
ducting behavior of samples of reduced dimensionality £ 
they can make electron pairing and long-range phase 
coherence occur at different temperatures in uncon- 
ventional superconductors^ lead to the Bcrczinskii- 
Kosterlitz-Thouless transition^ in two dimensions (2D), 
cause the destruction of long range phase order in in- 
finitely long one-dimensional (ID) wires^ and determine 
the resistive properties of ID wires of finite lcngthi 7 - - — 

We study the properties of superconducting ID wires 
in a model system: uniform isolated aluminum rings. An 
applied magnetic field threading a mesoscopic supercon- 
ducting ring drives transitions between states with differ- 
ent phase winding number (or number of fluxoids) . The 
process of changing the fluxoid number by 2ir is called a 
phase slip. Changes in fluxoid number also result in tran- 
sitions between states with different angular momentum 
(or vorticity) with results in jumps in measurable quan- 
tities such as the current around the ring. The nature of 
superconducting fluctuations in rings has generated sig- 
nificant interest. Fluxoid dynamics in individual rings 
have been probed as a function of ring size^^— mag- 
netic field ) 12 ' 14 ! 15 and temperature! 13 ' 16 ' 17 The occupa- 
tion of metastable fluxoid states has also been measured 
to determine a crossover from ID to 2D behavior in wide 
ringsi 18 ' 19 Phase slip rates have been studied in both con- 
ventional low T c [2(| and unconventional high T c [2l| su- 



perconducting rings. Ring inhomogeneities, such as weak 
links or nonuniform widths, have been studied as phase 
slip sites that can impact the ring's current-phase rela- 
tionship and fluxoid transitions! 15 ' 22 ~ — 

Transport measurements have long been used as a 
probe of superconducting fluctuations! 27 ' 28 These exper- 
iments measure a voltage that is directly related to the 
rate at which phase slips occur. In contrast, magnetic 
measurements such as ours are sensitive to the thermo- 
dynamic equilibrium current in the ring. Rather than 
tracking individual fluctuation events, this paper focuses 
on the effects of superconducting fluctuations on the 
ring's equilibrium supercurrent, /, as a function of ap- 
plied flux, $ a , measured in a temperature range near the 
critical temperature, T c . In order to reach equilibrium, 
phase slips must occur much faster than the measurement 
time. We focus on a regime, determined by the ring's 
physical parameters, where the distribution and switch- 
ing between fluxoid states dominate the ring's equilib- 
rium response to an applied field. Fluxoid states are the 
metastable superconducting solutions which correspond 
to local minima of the Ginzburg-Landau free energy func- 
tional and can be labeled by an integer phase winding 
number. In thermal equilibrium, when the ring's limited 
thermal energy is sufficient to cross the barrier between 
fluxoid states but still small compared to the saddle point 
energy, it spends almost no time near the saddle point 
solutions that allow phase slips. As a result, our equilib- 
rium measurement of the ring current depends only on 
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the energy difference between ffuxoid states. 

Theoreticai work using Ginzburg-Landau (GL) theory 
has predicted the current in the presence of an applied 
flux threading the ring in the opposite limit, where fluc- 
tuations between states with different ffuxoid numbers 
are inadequate to describe the response in the presence 
of all possible fluctuations including local phase and am- 
plitude fluctuations. Ambegaokar and Eckern applied a 
Gaussian approximation to GL to predict a mesoscopic 
persistent current driven by superconducting fluctuations 
above T c ) 29 i 30 which decreased exponentially with in- 
creasing T on the scale of the correlation energy. How- 
ever, the Gaussian approximation, accurate far above T c 
where the quadratic term in the GL free energy dom- 
inates, diverges as T approaches T c . von Oppen and 
Riedel used a transfer matrix approach to GL theory that 
accounts for all thermal fluctuations above and below T c 
to calculate the supercurrent and correct the divergence 
at T c £ Schwiete and Oreg then proposed a simplifica- 
tion of the full von Oppen and Riedel (VOR) theory that 
makes an analytic prediction for the ring's susceptibility, 
dl/d<fr a , in the limit where the superconducting coher- 
ence length is on the order of the radius* 3 ^ which provides 
a simple alternative to solving the VOR model numeri- 
cally. 

Direct measurements of the ring current as a function 
of applied flux are useful because it provides access to 
the thermodynamic free energy through the derivative 
/ = —dF/d$ a . While there are also interesting features 
in the full flux dependence; 3 - 2 - the fluctuation response 
is well captured by the zero field susceptibility. Conse- 
quently, in this paper we measure the ring's zero field sus- 
ceptibility as a function of temperature, d/(r)/d<^|^ = o, 
where 4> = $ a /&o and $o = h/2e is the superconducting 
flux quantum. 

A number of different experiments have used suscep- 
tibility measurements to study fluctuations in individual 
superconducting rings. 13 ' 20 ' 23 ' 32 ' 33 Zhang and Price stud- 
ied the phase slip rate and susceptibility as a function of 
temperature in a single Al ring^S The ring's geometry 
and long mean free path favored amplitude fluctuations 
that were expected to support a susceptibility response 
above T c . However, the observed susceptibility signal was 
an order of magnitude larger than predicted by GL the- 
ory. Koshnick et al£^ measured the susceptibility of 15 
individual rings with long mean free paths as a function 
of $ a . All rings showed a fluctuation induced susceptibil- 
ity response above T c . The magnitude of the fluctuation 
signal was large in the Little-Parks region, where an ap- 
plied flux suppresses the ring's T c (<I> a ) below its zero field 
value, T c ($ a = 0). The susceptibility responses of all the 
rings in that measurement were well described by von 
Oppen and Riedel's complete GL models 

Dirtier rings or rings with narrower geometries should 
exhibit ffuxoid fluctuations. Instead of generating an en- 
hancement in the susceptibility above T c , ffuxoid fluc- 
tuations can suppress the superconducting response well 
below T c . This paper focuses on rings that are likely to 



experience ffuxoid fluctuations. We start by describing 
the different thermal fluctuations, namely phase slips and 
ffuxoid fluctuations, experienced by our rings in Sec. Ill Al 
We then establish the physical conditions that support 
ffuxoid fluctuations in Sec. Ill Bl We outline a theory-^ 
where a thermal distribution of ffuxoid states suppresses 
the rings diamagnctism (Sec. Ill C[) . Our theoretical anal- 
ysis concludes by comparing our ffuxoid only model to 
a complete theory that includes all thermal fluctuations 
in the GL framework 2 - (Sec. Ill Dp . where we find good 
agreement in rings where fluxoid fluctuations dominate 
the response. We finally discuss our measurement tech- 
nique in Sec. IIIII and present data from two sets of ring 
samples with different mean free paths, which we com- 
pare to the theoretical models (Sec. MI C|) . 

II. FLUCTUATION THEORY 

We begin by describing the different types of ther- 
mal fluctuations encountered in the GL formalism ap- 
plied to a ring geometry, paying particular attention to 
the difference between phase slips (processes that change 
the fluxoid number by moving through a saddle point 
in configuration space) and fluxoid fluctuations (fluctu- 
ations in fluxoid number based on the energy difference 
between fluxoid states). We then discuss the onset of 
thermal equilibrium by considering the phase slip theory 
of Langer and Ambegaokar— as formulated for ring s 20 ' 34 
in the context of our actual measurement. Finally, we in- 
troduce fluctuation models that illuminate how different 
types of fluctuations influence the ring's susceptibility re- 
sponse. These model include a mean field model with no 
fluctuations, a model that includes only fluxoid states, 
the VOR model including all thermal fluctuations in the 
GL framework, and a harmonic oscillator approximation 
of the VOR model that includes only quadratic fluctua- 
tions. By comparing these models we can determine the 
type of fluctuations that dominate the response at dif- 
ferent temperatures and for rings with different physical 
parameters. 



A. Types of Fluctuations 

In the context of GL formalism we introduce a complex 
order parameter, i^(t), with an associated amplitude and 
phase. Fluctuations in the amplitude and phase of the 
order parameter are deviations in ip from the mean field 
solutions corresponding to local minima of the GL free 
energy. Fluctuations play a large role when the thermal 
energy of the system allows multiple wavefunctions to 
contribute to the ring's properties. 

The ring geometry of our samples imposes a con- 
straint on the order parameter phase. The periodic 
boundary condition requires that the phase be single val- 
ued, meaning it must wind by an integer multiple of 2tt 
around the ring. A phase slip is the process of changing 
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the phase winding number by 2tt. They usually occur 
by briefly suppressing superconductivity in a coherence- 
length-sized section of the ring.— Rings whose circumfer- 
ence, L, is longer than the superconducting coherence 
length, £(T), have multiple metastable states that differ 
by a phase winding or fluxoid number and their homoge- 
neous superfluid density. In a temperature range where 
phase slips occur within the measurement time, the ring 
fluctuates between its minimum energy fluxoid state and 
the metastable fluxoid states.— Such occupation of mul- 
tiple fluxoid states or fluxoid fluctuations can contribute 
to a suppression of the ring's superconducting response. 



B. Phase Slips and Equilibrium 

The ring's measured response represents its thermal 
equilibrium response when phase slips occur at a rate 
that is fast compared to experimental time scales (~ 
10ms), and no hysteresis is observable. In this section 
we derive a condition for the onset of phase slips and a 
separate condition for the onset of fluxoid fluctuations. 

We study switching between states driven by phase 
slips by adapting the theory of Langer and Ambegaokarfi 
Langcr and Ambcgaokar's treatment assumes a ID su- 
perconductor where the cross-section is smaller than the 
superconductor's coherence length, £(T), and penetra- 
tion depth, A(T), so ip has no radial variation. This as- 
sumption is particularly useful since exact solutions that 
minimize the GL free energy exist for ID systems. Addi- 
tionally Langcr and Ambcgaokar consider a homogeneous 
superconductor making phase-slips equally likely at any 
point along the ring's circumference. Langer and Ambc- 
gaokar's theory with a correction to the prcfactor by Mc- 
Cumber and Halperins predicts a phase slip rate by cal- 
culating the lowest energy pathway between two fluxoid 
states as defined by the energy barrier for the saddle point 
in wave function configuration space. Each fluxoid state, 
with energy -F m i n and phase winding 2im, represents a 
stable local minima of the GL free energy functional. The 
saddle point energies, -FLaddie, being stationary points of 
the free energy, must also satisfy the GL equations; how- 
ever, these solutions represent unstable configurations. 
The phase slip rate, T, then depends exponentially on 
the energy barrier AF(T, 0) = F saddlc (T, <f>) - F min (T, 0). 
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The prefactor, CI = ( J L/C)(AF/fc B T) 1 / 2 /r, was corrected 
by McCumbcr and Halpcrin based on arguments us- 
ing the time dependent Ginzburg Landau model^ where 
t = 7rft,/8fcs(T c — T) is the relaxation time. For phase 
slips to occur at a rate of 100 Hz, a rate that is approx- 
imately equal to our measurement time, the energy dif- 
ference AF/UbT may not exceed ~ 20. Although the 
prefactor is large, fl w 5 x 10 11 — 5 x 10 12 Hz, the expo- 
nential dependence overwhelms changes in the prefactor 



which can therefore be ignored. We find an approximate 
condition for the onset of phase slips from calculations of 
the energy barrier AF. 

We find expressions for -F m i n and -F sa ddie which repre- 
sent stationary points of the GL free energy. One dimen- 
sional Ginzburg-Landau theory introduces the GL free 
energy functional in the presence of a magnetic field rep- 
resented by the vector potential A as 
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a and (3 both depend on T, and a 2 / (3 = B c (T) 2 /fi is 
related to the superconducting critical field, B C (T). e* 
and to* are the charge and mass of the Cooper pairs and 
Ho is the permeability of free space. 

We look for stable solutions that locally minimize the 
Ginzburg Landau free energy functional. In a homoge- 
neous one dimensional ring fluxoid states have free ener- 
gies 
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where the critical field, and the ring volume (V = 
2i:Rwd) determine the ring's total condensation energy 
(F C (T) = VB c (T) 2 /2no). w is the ring width and d is 
the thickness. The dependence on £(T)/i? accounts for 
the suppression of the superfluid density by the phase 
gradient around the ring with coherence length £(T) and 
radius R. The Aharonov-Bohm flux, <fi = & a /&o, can be 
transformed into a shift in the boundary conditions for 
a wave function in a ring^ and therefore contributes to 
the energy in the same way as n. 

To find the saddle point energies, Zhang applied 
Langer and Ambcgaokar's phase slip theory to a ring 
geometry^ and by used the approximation L ^> £ (T) to 
calculate the saddle point energy 
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where L = 2ixR is the ring's circumference and S(T, (f>, n) 
is the normalized difference between the square of the 
order parameter amplitudes near and far from a phase 
slip event. 6(T, <f>, n) is a real number between and 1 
that satisfies the relation 
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For (b = n + 1/2, 5 = 1. 
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We are interested in a regime where L ^> £(F), and 8 
remains close to one for moderate n. Using the substitu- 
tion k = yl — 8 and expanding to lowest order in k, we 
arrive at a simplified expression for 8. 



8{T,<t>,n) = 1 



' y/3ir(2n - 2</> - 1)" 



L 



2^2 
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Returning to the energy expression, Eq. (|4|), we find 
minima of F sa ddie at n— 0=1/2 



min(F sad dio) = F C (T) 



'SV2C(T) 



I 



- 1 



(7) 



The energy expressions for F min (Eq. {3J) and F sa ddie 
(Eq. (j4j) are plotted as a function of flux for differ- 
ent fluxoid number n in Fig. Q] The minimum value of 
^saddle given in Eq. ([7]) is also plotted. The plots demon- 
strate how an increase in the ratio F/£(F) decreases the 
variation of both F mnl and F sa ddie with applied flux. For 
L/£(T) = 30 the variation in AF(</>) between </> = and 
4> = 1/2 is 30% of AF((j) = 0). This variation decreases 
to 12% for L/£(F) = 80. 
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FIG. 1: (Color Online) Free energy of fluxoid states (blue 
solid line, Eq. [3| and saddle point energies (black dashed 
line, Eq. [4]) as a function of applied flux, in units of the con- 
densation energy, F c . The red dotted line is the minimum 
value of Fiddle which is approximately equal to 3.8 times 
the condensation energy for a ring section of length f. (a) 
The respective energies when L/£(T) = 30. (b) The respec- 
tive energies when L/£(T) — 80. The energy barrier between 



adjacent fluxoid states, AF* 1 ' = F™^* 1 - 
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ferently than the saddle point energy AF. 



Fmin , scales dif- 



As L/£(T) becomes increasingly large, in rings 
with long circumferences, the flux dependence of both 
F sa ddie {T, (j)) and F m i n (T, 4>) becomes less important to 
the energy barrier AF(T, </>). This decreasing depen- 
dence on flux justifies the approximation that nei- 
ther F m j n nor F sa ddie depends on the flux in rings 
where L ^> £(F). Approximating F sa ddie(F, 4>) m 
min(Fsaddie) and F min (T, (f) « min(F min ), leaves 
AF(F) w 3.8wd£,B c (T) 2 /2fi , which is also independent 
of flux. The exponential dependence of the phase slip rate 
combined with the large attempt frequency in Eq. ([T]) 
imply phase slips occur when AF(T) < 20fceF. Qual- 
itatively this result makes sense because it states that 
phase slips occur when there is sufficient thermal energy 
to suppress the order parameter on the scale of the co- 
herence length. This energy scale differs from the ring's 
total condensation energy, F c , which is proportional to 
the volume of the entire ring. 

In contrast to the phase slip energy we just derived, the 
we obtain the energy associated with each fluxoid state 
by expanding the mean field GL free energy expression, 
Eq. ©, to lowest order in i{T)/R. 
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where 



I°(T) 
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%H R 2 
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(9) 



We introduce the parameter 7(F) which describes the 
inverse phase stiffness of the ring. Rings with a larger 
7(F) are more likely to exhibit fluxoid fluctuations. 7(F) 
is discussed in more detail in section III Dl 

When the energy difference between the lowest energy 
fluxoid states, AF^° = F^ 1 - F™=°, is comparable 
to the temperature, multiple fluxoid states contribute to 
the ring's equilibrium response. 

The criterion for phase slips (AF(F) w 
3.8wd£B c (T) 2 /2/j.q) scales differently with ring ra- 
dius compared to the energy difference between fluxoid 
states (AF^f(T) w ATtwd£B c {T) 2 /2/x (£/ ' R)) that de- 
termines onset temperature for fluxoid fluctuations. As 
stated above phase slips begin when AF(F) < 20fcsF, 
while fluxoid fluctuations onset for AF < 6fcsF. Sec 



section III CI Eq. (TT3"| . These two types of fluctuations 
onset at the same temperature when £(T)/L = 0.014. 
Therefore, in the limit where L 3> £(T) (specifically 
when £(F)/L < 0.014), the condition for fluxoid fluc- 
tuations is already satisfied when phase slips begin to 
occur. By applying the fluxoid condition, Eq. (|13p . and 
making an approximation that F « T c we can remove 
the temperature dependence and rewrite the condition 
as 7(F C ) > 16, 000. As a result, in the longest wires with 
the largest 7(F) a suppression of the response due to 
fluxoid fluctuations is expected as soon as the ring is 
able to reach thermal equilibrium. 

In the opposite limit, for small clean rings, if the con- 
dition AF^' < 6k B T is not met by the time AF(F) tn 
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k B T then thermal occupation of the saddle point solu- 
tions and other higher energy states become significant. 
This occurs in rings where 7(T C ) < 40. When fluxoid 
fluctuations onset they will not dominate the fluctuation 
response and their effect on the zero-field susceptibility 
is negligible. 

Using Eq. Q] we can predict the onset of phase slips in 
aluminum rings. A ring with a mean free path t e = 4nm 
has A(0) « 800 nm and £(0) w 85 nm5 Using standard 
temperature dependencies for these parameters, £(T) oc 
(1 - t)- 1 / 2 and \(T) oc (1 - i)" 1/2 where t = T C /T^L a 
ring of this material with w = 80 nm and d = 40 nm will 
have switching occur on experimental time scales down 
to 1.1 K when T c w 1.24 K, while a ring with d = lnm 
will continue to have phase slips down to 0.45 K. 

We have shown that the onset temperature for phase 
slips, AF oc £(T), and for fluxoid fluctuations, AF^ oc 
£(T) 2 /R, scale differently with the ring radius. How- 
ever, the presence of phase slips is a precondition for 
fluxoid fluctuations. The phase slip rate, Eq. (JTJ, de- 
pends on the temperature and the energy difference be- 
tween the saddle point energy and the minimized GL 
energy. In the L ^> £(T) limit, the rate is indepen- 
dent of applied flux and the ring circumference, F oc 
cxp(-3.8wd£(T)B c (T) 2 /2/.i Q k B T). This phase slip rate 
and the measurement time sets our criterion for equilib- 
rium. Namely a ring's measured response represents its 
equilibrium response when phase slips occur faster than 
the time over which the flux is swept through the transi- 
tion region. 

In the following sections we explore the effects of fluc- 
tuations on the ring's response in thermal equilibrium. 
Each of the models presented below includes a different 
set of fluctuations. By comparing the model predictions 
for different ring parameters we can pinpoint the effect 
of different fluctuations on ring response and set a phys- 
ical regime where each of the fluctuations will dominate. 
Specifically we find that in rings with weak phase stiffness 
a model including only fluxoid fluctuations accurately re- 
produces the ring response except for an apparent shift 
in T c indicating that fluxoid fluctuations dominate. 



C. Fluxoid Number Distribution Model 

We start with a model that includes only fluxoid fluc- 
tuations: fluctuations caused by integer winding or un- 
winding of the phase of the order parameter. This model 
is not complete because it does not include nonuniform 
phase fluctuations or amplitude fluctuations. Put an- 
other way this model includes only the large fluctuations 
between local minima in the GL free energy (see Fig. [TJ, 
and ignores all the small fluctuations about each local 
minimum as well as the saddle points and intermediate 
states. It is instructive to develop this model because 
comparisons between this fluxoid only model and more 
complete models shed light on what portion of the fluc- 
tuation response of a ring is due to solely to fluxoid fluc- 



tuations. As discussed in the previous section, we expect 
this model to be valid for rings with j(T c ) > 40. 

We return to the mean field ID GL free energy, Eq. ((3|), 
which is related to the ring current by I = —dF/d$ a . 
Taking the derivative yields an expression for the ring 
current of a state with n fluxoids. 



In(T, 



I°(T)(cf>-n) 1 



R 2 



(<f>-n)' 



(10) 



where I°{T) was given in Eq. ©. I a {T)(cf> - n) is the 
Mcissncr response which decreases linearly with increas- 
ing temperature close to T c , while the cubic term arises 
from pair-breaking. 

The energy associated with each fluxoid current state, 
F™ in , was derived in the previous section, Eq. flSJ. If 
phase slips occur at a high enough rate, so that the 
metastable fluxoid states are in thermal equilibrium as 
discussed in the previous section, we can model^^ the 
resulting current response as a Boltzmann distribution of 
fluxoid states. 



If(T,< 



E n In(T, <j>) exp (-FZJP, cj>)/k B T) 
E„exp(- J F« in (T,0)/fc s T) 



(11) 



We call the total ring current generated by fluxoid states 
If to distinguish it from the total ring current including 
all fluctuation states that will be presented in the next 
section. The susceptibility response of the ring at zero 
applied flux, dI(T)/d<fi\$ = o gives us access to the ther- 
modynamic free energy. In our rings £(T)/L < 1, so 
we make the approximation I n (T, 4>)fzI°(T)(<j> — n). The 
derivative of the total ring current at <f> = is 



dI F (T) 



= r(T) 1 - 



<t>=o 



J2 n 2t™ 2 exp(-crn 2 ) 
£„exp(-crn 2 ) 



(12) 



where a = I°(T)$ /2k B T. Eq. (TJ) shows that includ- 
ing a distribution of fluxoid states reduces the ring's sus- 
ceptibility response from the mean field value, 7°(T). 
However, the ring's supcrfluid density is not reduced. 
The second term in Eq. (fT2"|) is proportional to the RMS 
fluctuation of the fluxoid number, n. The magnitude of 
the reduction in susceptibility depends on a. When a is 
large, terms with n^0 are small and the susceptibility is 
approximately equal to the mean field value. When a is 
small, the n = ±1 terms begin to play a significant role. 
We define a criterion! for fluxoid fluctuations to reduce 
the Mcissner response by more than 5% when 



dI F{T) _ r0 ,^ , 12k BT 



r(T) < 



(13) 



as long as phase slips occur at a sufficiently high rate. 
We used this criterion in section III Bl to compare the on- 
set of fluxoid fluctuations and phase slips. In plots of 
the susceptibility vs. temperature we observed a sup- 
pression below the mean field value for susceptibilities 
below this cutoff. This downturn in the susceptibility 
signal, which occurs at T less than T c , is a hallmark of 
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the suppression of the diamagnetic response by fiuxoid 
fluctuations. More specifically, as stated in section III Bl 
we expect a visible downturn in rings with 7(T C ) val- 
ues between 40 and 16000. Rings with y(T c ) < 40 will 
not show a downturn as higher energy fluctuations over- 
whelm fiuxoid fluctuations. In rings with 7(T C ) > 16, 000 
fluxiod fluctuations will occur as soon as thermal equi- 
librium is reached. These rings will exhibit a response 
that is already suppressed below the mean field value for 
all data taken in thermal equilibrium. The next sections 
will be devoted to developing more complex models which 
include a more complete set of thermal fluctuations. 

Thus far, we have considered a fiuxoid model that pre- 
dicts the existence of the downturn in susceptibility be- 
low T c . In some rings, near T = T c , the R > £(T) 
assumption we made to obtain Eq. (jSJ) and (|12| is not 
strictly valid. As a result, the energy between succes- 
sive metastable states can no longer be approximated 
by the expansion in Eq. ([5]). By including the quar- 
tic term of Eq. the GL free energy vanishes rather 
than increasing indefinitely for 4> — n > L/£(T). Thus, 
the Boltzmann distribution, Eq. is also not well 

defined because summing over all n leads to a divergent 
denominator. The numerator on the other hand remains 
finite since states with </> — n > L/£(T) do not contribute. 
Furthermore, our treatment thus far has ignored phase 
fluctuations that are not uniform around the ring and all 
amplitude fluctuations. 



D. von Oppen and Riedel Model 

To address these issues, we compare our simple fiux- 
oid model to the model of von Oppen and Riedel^ VOR, 
which generates numeric solutions that include all ther- 
mal fluctuations within the GL framework in homoge- 
neous rings. Applying a harmonic oscillator approxima- 
tion to the VOR model, discussed in the next section, 
provides a direct mathematical connection between the 
VOR and fiuxoid model discussed in the previous section. 

Following von Oppen and Riedel^ we begin with the 
expression for the GL energy functional given in Sec. [TT] 
Eq. ([2]) . We map the free energy onto a one dimensional 
ring geometry with no lateral variation of the order pa- 
rameter so ip(r, z, 9) = ip(9) and dx 3 = wdRdO. We can 
then redefine Rd6 as dx. 

We rewrite Eq. (|2|) using reduced variables tp(x) = 
ip(x)^\a\//3, V = fV, and x = x£. f(T) is the su- 
perconducting coherence length and is given by £(T) = 
hj \J1m*u. 



F$(x)] = 
E (T)k B T 



A(T)/2 

'-A(T)/2 
27TI 



v\m\ 2 + ^\m\ 4 



A(T) 



4> ip(x) 



dx 



(14) 



r/ is +1( — 1) for temperatures above (below) the super- 
conducting critical temperature T c . A(T) is the reduced 
circumference A(T) = L/£(T) = y/8irk B \T - T c \/E c and 
E (T)k B T = wd4(T)B c (T) 2 /n is the condensation en- 
ergy of a ring section of length £(T). The correlation 
energy for the ring, E c = ir 2 hvf£ e /3L 2 , includes the 
mean free path, f, . and the fcrmi velocity, Vf, which is 
2.03 x 10 6 m/s in aluminum. Eq(T) can also be written 
as 



E (T) 



(2tt) 5 / 2 
21C(3) 



k B \T-T c 

E r 



3/2 



E r M, 



off 



k R T 



(15) 



where £(3) = 1.021 is the Riemann zeta function. M = 
k 2 wd/4ir is the number of transverse channels, kf is 
the fermi wave vector, which for an aluminum ring is 
kf = 1.75 x 10 10 m~ 1 . Including disorder results in an 
effective number of channels, M e g = M£ e /L. 

We obtain the thermodynamic expression of the cur- 
rent from the flux derivative of the ring's partition func- 
tion. 



I(T, 0) = -k B T 



_1 d_ 



(16) 



The partition function is the path integral of the GL 
energy. 



Z S c = I [dip{x)} [dip* {x)\ exp 



-F[m\ 

k R T 



(17) 



,.'SS 



to 



The VOR model uses a transfer matrix techniqu 
map the Ginzburg-Landau path integral partition func- 
tion, Eq. (|17|) . onto another partition function 

oo oo 

Z= exp(-*27rZ0) ^ exp(-2 J B (T)A(T)# n , i ) 

l — oo n—0 

(18) 

where S n> i are the eigenvalues of the fictitious 2D single- 
particle Hamiltonian, 



H = - 



1 



1 . 1 



8E (T) 2 2 ' 4 



(19) 



We define p = (2E (T)) 1 ^ 3 r and rewrite Eqs. (15]) and 
(UHJ) to emphasize the parameter j(T)M- 



Z= cxp(-i27r/0)^exp(-7(T) 1 / 3 ( l' lM )(2O) 

(21) 



n=0 



H 



1 a , 1 A(T) 2 „ , 1 . 



2 v ' 2 7 (r) 2 /3 P A 9 



The temperature dependence is set by the coherence 
length through the relation A(T) = L/£(T). The pa- 
rameter 



7C0 



A(T) 3 42C(3) k B T 



2E (T) 7T M eS E c 



(22) 
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can be thought of as a measure of inverse phase stiffness, 
which determines the type of fluctuations that dominate 
the ring's susceptibility response.— The definition of 7 
introduced in Koshnick et al££ made the approximation 
T w T c . The larger temperature range explored in this 
paper makes it necessary to reintroduce the T depen- 
dence. We use the relation for J°(T) given in the second 
part of Eq. © to compare the VOR model to the mean 
field and fluxoid models. 

Eqs. ([20]) and pTj) can be solved numerically. The 
Hamiltonian can be rewritten as a harmonic oscillator 
with a quartic perturbation. We can then write matrix 
elements in terms of the coefficients and diagonalized nu- 
merically to find the eigenvalues 1 34 ' 39 The eigenvalues are 
used in the partition function, Eq. (|20|) . and substituted 
into the thermodynamic equation for the current, Eq. 
([T5|) . to generate the full current response. We find the 
zero-field susceptibility by taking a derivative with re- 
spect to applied flux at <f> = 0. 

Analytic solutions can be instructive, and as a result 
it is useful to find approximations to the full VOR model 
that are valid over some set of ring parameters or temper- 
atures. One such approximation is to ignore the quartic 
perturbation to the Hamiltonian, which then takes the 
form of a simple harmonic oscillator. 



E. Harmonic Oscillator Model 

The harmonic oscillator (HO) approximation is valid 
at temperatures well below T c , where the wave functions 
contributing to Eq. (|20j) only extend over a narrow re- 
gion around the minimum of the Mexican hat potential 
of Eq. pTj) . so that the latter can be approximated by 
a quadratic expansion. In this case, fluctuations from 
the quartic nature of the potential should not play a sig- 
nificant role. We refer to the fluctuations in this model 
as quadratic fluctuations, rather than Gaussian fluctua- 
tions, to avoid confusion with small order parameter fluc- 
tuations above T c , which arc often referred to as Gaussian 
fluctuations. 

Eigenstates have the form r = \r\ exp(zZ</>), so Eq. J2~T|) 



(a) 



can be written as a ID problem, H = —\-§pz + V(r) 
where 



I 2 1 A(T) 2 



2r 2 2 7 (T) 2 /3 



— r 

4 



(23) 



Expanding V(r) about its minimum at R m (l) leads to 
the eigenvalues 
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FIG. 2: (Color Online) Theoretical susceptibilities calculated 
using the mean field model, black dotted line Eq. ((9}, the flux- 
oid model, purple dot-dash line Eq. (|12|l , the VOR model, red 
solid line Eqs. (|16ll20ll2"Tj) . and its approximate HO solution, 
green dashed line Eqs. (|16II26|I . for rings with w — 80 nm and 
T c = 1.55 K. a) 7(T = T c ) = 3. The VOR model predicts a 
susceptibility above T c . b) j(T = T c ) = 3100. A downturn 
occurs at dl/d(j> « 12fc B T/$o « 120nA, T w 1.52,K. The 
fluxoid, HO and VOR models reproduce the overall lineshape 
of the downturn, up to an offset in T c . However, the three 
models predict downturns of different sizes with the largest 
predicted by the VOR model, c) j(T = T c ) = 170, 000 Flux- 
oid fluctuations dominate the response over a wide tempera- 
ture range and the HO and fluxoid models become increasing 
accurate predictors of the full fluctuation theory. For all val- 
ues of j(T = T c ), the VOR and HO response well below T c 
only match the mean field and fluxoid predictions if T c is 
renormalized. 



I 2 



2R m (iy 



i?m(0 4 



+ w(n+l/2) (24) 



where w = ^A^fJWW + 3i?m(0 2 + WjRJff. 

Only terms that change with I, the angular momen- 
tum coordinate in the fictitious Hamiltonian, contribute 
to the flux dependence of the partition function, thus 



only these terms contribute to the thermodynamic ring- 
current. If we make an approximation and only include 
the Z 2 /2i? m (0) 2 terms, where i? m (0) is the value for r 
that minimizes V(r) when I = 0, the current from Eq. 
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QB is 



I ho (T, i 



shl (2tt2</>) exp(/ 2 7 (T) /2A(T) 2 ) 

fcs? 1 ;=i 



$ 



1 + ]T 2 cos {2nl<t>) exp(7 2 7 (T)/2A(T) 2 ) 



,(25) 



which is exactly equivalent the fluxoid current shown in 
Eq. Through this approximation we are able to 

show a direct link between the harmonic oscillator ap- 
proximation to the VOR model and the fluxoid model. 
Including the second two terms of Eq. (|24|) . which ac- 
count for the angular momentum dependence of lu and 
R m (l), we get 



Z= J2 cxp(-z2^0)exp(- 7 (T) 1 /V(i? m (O)) 

I — — OO 

cxp(- 7 (T) 1 / 3 a;/2) 



1 -exp^^T) 1 ^)' 



(26) 



Using this simplified partition function we can find the 
ring's current and consequently its susceptibility in the 
limit where we ignore only quartic fluctuations. 



F. Comparison of Models 

We have presented the theoretical basis for four models 
including: the mean field model, the fluxoid model, the 
harmonic oscillator model and the von Oppen and Riedel 
model. We now compare the physics captured by each 
model by plotting the theoretical susceptibility response 
predicted by each model as a function of temperature for 
rings with three different 7 (T = T c ) parameters in Fig. 

El 

The mean field model is our baseline. It gives the ring 
response in the absence of all superconducting fluctua- 
tions. At the other extreme, the VOR model incorporates 
all thermally activated superconducting fluctuations into 
its derivation of the ring response. In between we have 
the fluxoid model, which includes only fluxoid fluctua- 
tions and the harmonic oscillator model. By comparing 
these models for rings with different 7 (T = T c ) we can 
get a sense of which fluctuations dominate the response. 

One striking feature in Fig. [5] in all three plots for all 
values of 7 (T = T c ) 1 is that both the VOR model and its 
HO approximation have an offset in the linear regime, far 
below T c , compared to the mean field or fluxoid model. 
This downshift in T c appears to reflect a renormaliza- 
tion in T c due to consideration of all possible fluctuation 
modes. 

Fig. [!» shows a ring with 7 (T = T c ) = 3. The low 
gamma parameter means it has a strong phase stiffness. 
The susceptibility at T = T c is 210 nA. We get this num- 
ber by looking at the value of the VOR model at T c in Fig. 



[2ja). We compare this number to 12/c B 77$ ~ 120 nA, 
for T = 1.5 K predicted as the point where we expect 
a downturn in susceptibility due to fluxoid fluctuations, 
Eq. ([TU]) . Since the fluxoid criterion is smaller that the 
susceptibility at T c , susceptibility enhancing amplitude 
fluctuations at and above T c overwhelm the susceptibility 
reduction expected from fluxoid fluctuations. A down- 
turn is not observable, instead the small 7 (T = T c ) leads 
to a susceptibility signal above T c . 

When 7 (T = Tc) = 3100, as shown in Fig. HJb), 
the fluxoid induced downturn becomes visible below T c 
starting at T m 1.52 K and 120 nA, as predicted by our 
fluxoid criterion, Eq. (|T3"|) . All three fluctuation mod- 
els qualitatively reproduce the shape of the susceptibility 
suppression. As expected, VOR model predicts a greater 
susceptibility suppression than the fluxoid or HO mod- 
els, because only the VOR model includes all thermal 
fluctuations. The excess suppression between the fluxoid 
and VOR is presumably due to contributions from non- 
homogeneous phase winding solutions, amplitude fluctu- 
ations, or both. While the excess suppression between 
the HO and VOR models is due to fluctuations caused 
by the quartic nature of the potential. 

For 7 (T = T c ) = 170, 000, shown in Fig. [^c), the sus- 
ceptibility response is dominated by fluxoid fluctuations, 
shown by the almost identical lineshapc shared by the 
fluxoid model and the VOR model. The total response is 
also well represented by the harmonic oscillator approxi- 
mation showing that in this region nearly all fluctuations 
are quadratic in nature. 

Fig. Htc) shows a larger temperature range than the 
previous panels, and the GL approximation that T is 
close to T c is not valid over the whole plot. GL theory 
is strictly valid in the range of temperature where the 
linear mean field response approximates a temperature 
dependence that goes as (1 — i 4 ), t = T/T c , shown as an 
orange dotted line. An alternative criterion is that T > 
A(T), where A(T) is the superconducting gap. These 
both result in approximately the same range of validity. 
GL theory has been applied with success at temperatures 
far from T c , but interpretation of results in this regime 
should be treated with caution. The (1 — t 4 ) dependence 
is not included in panels (a) and (b) because all plotted 
temperatures lie within the valid range. 

In the next section we describe our measurement of 
ring susceptibility for rings with different 7 (T = T c ) val- 
ues. Wc find good agreement between our data and the 
fluctuation response predicted by the fluxoid and VOR 
models. 



III. SAMPLE AND MEASUREMENT 
TECHNIQUE 

A. Sample Preparation 

We measured quasi-one-dimcnsional superconducting 
rings in a dilution refrigerator— with a scanning SQUID 
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susceptometer— that was explicitly designed for this pur- 
pose. We focus on data from two different samples ex- 
pected to exhibit fluxoid fluctuations. Sample I's rings 
were fabricated and measured previously^ The rings 
were narrow and dirty with T cl m 1.5 K. They were made 
by depositing a 40 nm thick Aluminum film by e-beam 
evaporation at a rate of about lA/s and a pressure of 
approximately 10 -6 mBar on a Si substrate patterned 
with poly(mcthyl methacrylate) (PMMA) resist. During 
the deposition, the rate temporarily dropped to a negli- 
gible level for about 10 min and subsequently recovered. 
This delay caused the formation of two superconducting 
layers separated by an AlOx tunneling barrier. The cou- 
pling between the two Al layers depended on the width 
of the rings with narrow rings (w < 190 nm) and wide 
rings (w > 250 nm) showing a single order parameter. In- 
termediate widths showed evidence of weak interactions 
between the two layers leading to two order parameter 
effects^ In this work we only present data from the nar- 
row rings which showed no two order parameter behavior. 
However, due to the oxidization process we suspect the 
thinnest rings have a large oxidized layer that reduces 
the thickness of the superconductor. Consequently we 
expect that these rings have an effective height that is 
less than 40 nm. We can test this prediction by extract- 
ing the ring's cross-section from fits to the VOR model. 

The rings on sample II were fabricated specifically for 
this paper. The fabrication process was almost identi- 
cal to the rings from sample I except the evaporated 
film was thinner, d = 15 nm, and there was no inter- 
ruption in the evaporation. The deposited rings were 
wide and dirty with T cU « 2.1 K. Of the many fabri- 
cated rings of different widths and radii, only the widest 
rings, w w 850 nm, had a diamagnetic response. The 
next widest rings, w « 450 nm, showed no signs of super- 
conductivity indicating that they were oxidized through- 
out. This evidence makes it difficult to predict with cer- 
tainty what portion, if any, of the 850 nm rings are also 
oxidized. Although the ID approximation, w >> £(T), 
is no longer strictly valid for these rings the theory still 
provides a good fit to our data. For each sample we used 
Ginzburg Landau modcl o 20 ' 33 i 37 to fit a zero temperature 
penetration depth Ai(0) sa 800 nm, An(0) ~ 1-5 /im, and 
coherence length £i(0) « 80 nm, £n(0) « 30 nm. 



B. Measurement 

Our measurements arc done with a voltage biased 
DC SQUID susceptometer amplified by a series-array 
SQUID preamplifier 42 The SQUID is mounted on a 
piezo-resistive scanning assembly^ which is connected to 
the mixing chamber of a dilution refrigerator through a 
single copper braid. The temperature of the scanner and 
sample is controlled with sub-milikelvin precision though 
feedback. The SQUID sensor's counter-wound geome- 
try, with on-sensor modulation coils for feedback, enable 
cancellation of an applied field to one part in 10 4 4i The 
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FIG. 3: (Color Online) Plots of SQUID response vs. applied 
flux at different temperatures for a ring from each of the two 
samples. Rings b and e refer to specific rings plotted in Fig. 
H Ring b T c = 1.56 K and ring e T c = 2.08 K. The curves 
evolve from non-hysteretic with no fluxoids at low tempera- 
tures, through a hysteretic regime, to non-hysteretic with a 
change in fluxoid number at every applied flux quantum near 
T c . The orange dotted lines are fits to the GL current, Eq. 
(|10[) at low temperatures and to the Boltzmann distribution, 
Eq. JTTJ , at high temperatures. We extract the ring's suscep- 
tibility at each temperature by taking the derivative at cj> = 0. 



ring current is measured by positioning the SQUID about 
1 ixm above the ring and recording the flux induced by the 
ring's current in the SQUID's 4.6 fjxn diameter pick-up 
loop. During the measurement, the applied flux thread- 
ing the ring is varied by several flux quanta at a few Hertz 
by an on-sensor field coil. This measurement is repeated 
13 /im above the ring and the ring signal is computed as 
the difference between the two positions for each value of 
applied flux. This procedure allows us to achieve an ad- 
ditional three orders of magnitude of background cancel- 
lation. A more detailed description of the measurement 
system was given by Koshnick et al£^ 

We plot the flux induced in the SQUID's pick-up loop 
as a function of the flux applied by the field coil in Fig. 
[3] for two different rings. The measurement is repeated 
to record the full temperature dependence of the ring's 
response. The ring current, /, is coupled as flux into the 
SQUID pick-up loop through the mutual inductance, M. 
'J'squid = MI. We estimate the mutual inductance be- 
tween the SQUID pick-up loop and a ring by calculating 
the mutual inductance between two on axis rings with 
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radii r% and a distance z apart. 

/>oo 

M = ir(J,orir 2 / cke~ K|z| Ji(/cn) Ji(«r 2 ) (27) 
Jo 

Ji is a Bessel function of the first kind. For all our 
mutual inductance calculations we assume a ring-pick- 
up loop separation of 1 //m. Through a separate fitting 
technique^ we estimate the actual distance between the 
pick-up loop and the ring to range from 0.75 — 1.1 //m. 
Ring currents and susceptibilities quoted later in this 
paper have error bars that reflect this systematic un- 
certainty in the coupling factor which would shift the 
entire dataset. For example, the ring-SQUID mutual 
inductance between a ring with a radius of 1.5 [im po- 
sitioned 1 /im above the SQUID's pick-up loop is M — 
0.718 $o/ n A. The error on the height from 0.75 — 1.1 /im 
gives M = 0.669 - 0.849 <V nA , but because the SQUID 
is kept a constant distance from the ring for all data 
points taken at the same time on the same ring this sys- 
tematic error in the height would shift all data points 
together. 

The ring response curves plotted in Fig. [3] evolve from 
cubic and non-hysteretic at low temperatures through 
a hysterctic regime to periodic and non-hysteretic near 
T c . At low temperatures the current response is well 
described by the GL current with no phase windings, Eq. 
(|10p with n = 0, shown as a orange dotted line in panels 
(a) and (e) of Fig. [3J As the temperature increases the 
applied flux causes the ring to transition to higher fluxoid 
states; however, the phase slip rate becomes comparable 
to the measurement time only at large applied flux, thus 
leading to a hysteretic response. 

Finally, as the temperature approaches T c the phase 
slip rate becomes fast compared to the measurement time 
at small applied flux and the ring relaxes to thermal equi- 
librium. Thermal fluctuations are strong enough for tran- 
sitions to occur within some range of 4> = to/2 where 
to is an odd integer. The ring's response is no longer 
hysteretic and can be approximately modeled as a Boltz- 
mann distribution of all fluxoid states, Eq. (fTT|) . shown 
as a orange dotted line in panels (d) and (h) of Fig. [3] 
We extract the magnetic susceptibility of the ring at each 
temperature by fitting a low order polynomial to obtain 
the slope at <fi = 0. 

C. Susceptibility Data 

We measured thirty-eight rings on sample I and twelve 
rings on sample II. Sample I was fabricated and mea- 
sured primarily for a different experiment^ As a result 
only eight of the rings measured have sufficient suscep- 
tibility data over a wide enough temperature range to 
test the theories presented in the previous section. Two 
representative rings were selected for this paper. The 
three rings from Sample II were chosen to show a vari- 
ety of ring parameters, and because they had the most 
dense susceptibility data over the important temperature 



range. The set of five rings allows us to explore the effects 
of ring size and cleanliness on the fluctuation response. 
Fig. 5] shows the susceptibility vs. temperature data 
for those five rings. Each of the ring's physical param- 
eters are given in table HI We extracted the ring radii 
from the flux periodicity of the ring's response in ther- 
mal equilibrium and confirmed the measurement though 
SEM imaging. The ring thicknesses were measured with 
an AFM, and the width with SEM. Fitting to the VOR 
model allowed us to estimate values for the ring's cross- 
section and mean free path. We used the measured ring 
width and thickness plus an additional error factor as an 
upper limit on the cross-section parameter in the VOR 
model for rings (c-e) . No lower limit was enforced due to 
the possibility of oxidation reducing the superconducting 
cross-section. 
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TABLE I: Table of ring values. Values for the cross-section 
and mean free path, extracted from fits to the VOR model, are 
given with their 95% confidence interval. An upper limit of 
14790 nm 2 was enforced on the ring cross-section to constrain 
the fits for rings (c-e). 
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TABLE II: Table of fitted values. We used the temperature 
independent portions of r y(T) and Eq(T) as fit parameters in 
the VOR model. This table reports values for j(T = T c ) and 
Eo(T = 0) as well as the limits of the 95% confidence interval 
obtained from bootstrap analysis for the data presented in 
Fig- HI 

Fig. El plots the susceptibility vs. temperature curves 
for five rings. The blue susceptibility data points repre- 
sent the slope at <f> — of the SQUID response at different 
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FIG. 4: (Color Online) Zero field susceptibility (blue dots) 
vs. temperature plotted for five different rings. In all plots 
the solid red line is a fit to the VOR model, Eqs. J20] [21)) . 
The dotted black line is the expected mean field susceptibility 
given by I (T), Eq. ([9|. The gray shaded area represents 
the 95% confidence interval from bootstrapping. The error 
bars represent the systematic uncertainty in the SQUID-ring 
mutual inductance. 



temperatures scaled by the ring-SQUID mutual induc- 
tance to get the ring current. The error bars represent 
height errors in our calculation of the mutual inductance, 
Eq. (f2~T)) . This error is systematic and expected to be the 
same for all points in a panel. Using T c , and the temper- 
ature independent portions of 7(T) and Eo(T) as the free 
parameters, the red line is a fit of the data to the VOR 
model, Eqs. ([201 EH]) . The fit results used to generate 
the red curves are given in table |TTJ We report values 
for 7(T) at T c and E (T) at T = 0. The reported T c 
represents the nominal mean field T c entering the VOR 
models The fitted values of j(T = T c ) are also listed on 
each of the plots. The black dotted line is the mean field 
ring response, Eq. (j9|), which one would expect if no fluc- 
tuations were present. Deviations in the data from the 
black dashed line show the influence of fluctuations on a 
given ring. Finally, the gray region of the curve is the 95% 
confidence interval (CI) obtained from bootstrapping. 

Using the fit results from tablc|TT]along with the known 
values of the ring radii given in table U we can extract 
values for the ring's cross-section and mean free path 
from expressions for E (T), Eq. (fT5|). and j(T), Eq. 
(|22[) . The ring parameters obtained in this way are give 
along with their 95% confidence intervals in table HI Due 
to the evaporation conditions discussed previously, we're 
not confident that the entire cross-section of each ring is 
superconducting. For the two rings on sample I the fit- 
ted cross-sections arc smaller than the values found using 
AFM/SEM, which confirms our suspicion that a portion 
of the ring is oxidized. The data from the three sample 
II rings is within the downturn region, ie the decrease in 
the susceptibility is not linear even at the lowest plotted 
temperatures. Practically we are limited on the low end 
of the temperature range by the point where the SQUID 
response curves go hysteretic. A three parameter fit is 
under-constrained and it is consequently difficult to get 
accurate VOR fits without susceptibility data at lower 
temperatures including the point where the data is re- 
duced from the linear response. 

As a result, for rings (c-e) we put a strict upper limit of 
wxd = 14790 nm on the cross-section, which acted as an 
additional constraint on the VOR fits. The cross-sections 
extracted from fits to the constrained VOR model for the 
rings on sample II agree well with the AFM/SEM cross- 
section indicating little oxidation. A similar limit was not 
applied to rings (a-b) because data in the linear suscep- 
tibility region kept the fit from being under-constrained. 

The agreement between the susceptibility data and fits 
to the VOR model are good for all rings except ring (d), 
where it is clear that the VOR model does not capture 
the shape of the data at high temperatures. It is unclear 
why the VOR model provides a poor fit for this ring. It 
is possible that errors from extracting the susceptibility 
near T c , errors that are not accounted for in the error 
bars, are particularly large for measurements on this ring. 

Looking at the sequence of five rings it is clear that the 
extent of the suppression of superconductivity increases 
as j(T = T c ) increases. This is just what we expect for 
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FIG. 5: (Color Online) Comparison of the four models plot- 
ted using parameters obtained from fitting the data from ring 
(c) to the VOR model. Fitting to the fluxoid model would 
yield an equally good fit with a slightly different T c . 



we made throughout this paper. Fluxoid fluctuations 
not only suppress the rings superconducting response 
but play an increasingly large role in the suppression 
as y(T = T c ) increases. We showed that our suscepti- 
bility vs. temperature data is well described by a GL 
model for homogeneous rings, developed by von Oppen 
and Ricdel^ that includes all thermally activated fluctu- 
ations, but the largest gamma rings can be equally well 
described by our simple fluxoid only model with a shifted 
T c . Furthermore we can use fits to the VOR model to re- 
produce some of the rings physical parameters including 
the cross-section and mean free path. Finally, by using 
VOR fit parameters we can employ our two approximate 
models, the fluxoid model, and the harmonic oscillator 
model, to determine the how much of the suppression is 
due to either fluxoid fluctuations or quartic fluctuations, 
done for ring (c) in Fig. 



a series of rings where fluxoid fluctuations play a larger 
and larger role. 

Ring (a) shows an enhancement of the superconducting 
response above T c . This response is caused by amplitude 
fluctuations and has been studied by Koshnick et al£^ 
and Zhang and Priced As we showed in our descrip- 
tion of the theoretical models only the VOR model can 
correctly reproduce the upturn in susceptibility above T c . 
The fluxoid fluctuation model cannot reproduce a current 
above T c and will fail for any ring with a 7(T C ) < 40. 

The remaining four rings in Fig. 0] show a suppres- 
sion of the susceptibility signal below the mean field 
response (shown in black dots). However, of the plot- 
ted rings only ring (b) has a large enough temperature 
range to observe a downturn from the linear regime. The 
full temperature range plotted for rings (c-e) is already 
deep in the suppression region. This is due to the fact 
that j(T c ) > 16,000 for rings (c-e). In the next section 
we expand the temperature range by adding susceptibil- 
ity data from lower temperature hysteretic ring response 
curves, to confirm that the response is suppressed from 
the mean field value. The downturn for ring (b) occurs 
at T w 1.52 K and 120 nA, which corresponds to the cri- 
terion for fluxoid fluctuations given in Eq. (|13[) . Such 
agreement validates our criterion for the onset of suscep- 
tibility suppression driven by fluxoid fluctuations. 

We have shown that the VOR model describes the tem- 
perature dependence of the susceptibility. To get a feeling 
for the type of fluctuations that play a role in the ring 
response we plot the fluxoid model and the HO model in 
addition to the VOR model and mean field model for ring 
(c) in Fig. [5] It is clear that fluxoid fluctuations cause 
the majority of the suppression, with quadratic fluctu- 
ations described by the HO model contributing to the 
renormalization of T c and quartic fluctuations described 
by the VOR model playing only a minor role. In fact the 
fluxoid model would fit the data equally well with just a 
shift in the T c . 

The dataset, taken as a whole, confirms the points 



D. Hysteretic Susceptibility Data 

For rings (c-e) in Fig. [4] we expect the onset of the 
downturn induced by fluxoid fluctuations in a temper- 
ature range where the SQUID response curves are hys- 
teretic, as shown in Fig. |3l This is due to the fact that in 
these thinnest, dirtiest rings L >> £(T) and fluxoid fluc- 
tuations are already energetically favorable at the tem- 
perature when phase slips begin to occur, as discussed in 
Sec. Ill Bl Fluxoid fluctuations are never energetically fa- 
vorable for ring (a) and they onset well after phase slips 
in ring (b). Phase slips onset at ~ 1.3 K while fluxoid 
fluctuations onset at ~ 1.51 K. 

To demonstrate that the data presented represents a 
real reduction in the ring response we examine the sus- 
ceptibility signals at lower temperatures that fall in the 
hysteretic regime. We evaluate susceptibility in the hys- 
teretic regime by taking the slope at zero current on the 
long continuous sides of the hysteretic curves. These sus- 
ceptibility data points are shown as green dots in Fig. [5] 

Fig. [5] shows susceptibility data in the hysteretic 
regime (green points) and reproduces the susceptibility 
data from the non-hystcretic regime (blue points) from 
rings (b-e) in Fig. H) Fig. |H] also shows the (1 — i 4 ) 
dependence, plotted as an orange dotted line, that sets 
the validity of our GL based models. In Panel (a) the 
susceptibility in the hysteretic regime has dropped below 
the mean field prediction and instead closely follows the 
(1 — t A ) curve. 

Notice that the hysteretic data points in rings (c-e), 
panels (b-d) , follow the mean field curve until a crossover 
point when they line up with the VOR model and the 
higher temperature susceptibility data. This provides ev- 
idence that the susceptibility measured from the SQUID 
response curves in thermal equilibrium is suppressed from 
the mean field value. The drop in susceptibility from the 
mean field value occurs when phase slips occur at a suffi- 
ciently high rate and multiple fluxoid states compete to 
suppress the response. 
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IV. CONCLUSIONS 

Superconducting phase slips in one dimensional rings 
and wires have been the subject of theoretical and ex- 
perimental interest for decades. While phase slips in ID 
structures determine the onset of resistance, the fiuxoid 
processes we described here cause the loss of another hall- 
mark of superconductivity, the ability to screen magnetic 
field. In this paper we have outlined four models that de- 
scribe the effects of superconducting fluctuations on the 
susceptibility response in rings. We have shown that ring 
responses for rings with various physical parameters can 
be characterized by a model by von Oppen and Ricdcl 
for uniform rings that includes all thermal fluctuations. 
However, by comparing the models we can determine the 
types of fluctuations that contribute to the response of 
a given ring. We specifically found that for rings with 
weak phase stiffness the ring response can be described 
using a fiuxoid only model, indicating that these types 



of fluctuations are the dominant cause of suppression of 
the susceptibility signal. One could imagine extending 
this ring system to a weakly connected grid, linking our 
results to the field of percolation superconductivity. Ad- 
ditionally, achievable experimental conditions allow fiux- 
oid fluctuations to occur at temperatures down to 50 mK, 
which could provide an experiment setup for examining 
the quantum mechanical behavior of a ID ring4^ 
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FIG. 6: (Color Online) Susceptibility data from rings (b), 
(c), (d) and (e). The green points represent the slopes of the 
hysteretic curves, which estimate the susceptibility in the hys- 
teretic regime. The error bars account for a systematic error 
in the coupling constant that would shift all points together. 
The blue points give the susceptibility of the ring's response 
in thermal equilibrium. The red solid line is a fit of the blue 
non-hysteretic data to the VOR model and the black dotted 
line is the mean field response. Also plotted is the (1 — f 4 ) 
temperature dependence which places approximate limits on 
the validity of GL. 



